Crop–livestock integration enhanced soil aggregate-associated carbon and nitrogen, and phospholipid fatty acid

Integrated crop–livestock (ICL) production enhances diversification and provides ecosystem benefits by improving nutrient cycling and energy efficiency, thus, increasing overall farm productivity. However, a detailed study is needed to understand the influence of crop diversification and grazing animals on soil aggregation and associated carbon (C) and nitrogen (N), and microbial properties, especially compared with a grazed native pasture. We investigated the soil aggregate size distribution and associated C and N fractions, glomalin-related soil protein, and soil phospholipid fatty acid (PLFA) to understand the collective influence of livestock grazing of crop residue and cover crops (CC) and compared it with native pasture and non-grazed traditional production systems. The study was conducted in South Dakota at four different locations consisting of three long-term (> 30 years) on-farm sites: 1 (Salem), 2 (Bristol), 3 (Bristol) with three treatments that included ICL (corn, Zea Mays L.-soybean, Glycine max L.-oats, Avena sativa L.-CC with cattle grazing); natural ecosystem (NE) or native pasture; and control (CNT) (corn–soybean-without grazing or CC). Experimental site 4 (Beresford) with study duration of 3-year consisted of oats, oats with CC, oats with CC + grazing, and grazed pasture mix. Soil samples were collected from 0 to 5 cm depth at all four sites in summer 2019. Data showed that at sites 1 and 2, ICL had significantly (P ≤ 0.5) greater fractionation of 0.053–0.25 mm and > 4 mm aggregates compared with NE and CNT. At site 1, ICL showed significantly higher soil organic carbon (SOC, 36–49%) and higher nitrogen (33–44%) in > 4 mm aggregates than NE and CNT. At site 2, ICL had 32–41% higher SOC than NE and CNT for 0.25–0.5 mm aggregates. At site 1, NE enhanced total phospholipid fatty acid (PLFA), total bacterial biomass, gram (+), gram (−) bacteria than CNT, however, it did not vary significantly than ICL. Grazed pasture mix at site 4 had higher total PLFA (40.81 nmol g−1 soil) than the other treatments. The principal components 1 and 2 accounted for 33% and 22% of the variation, respectively, where the majority of the microbial compositions and aggregate-associated C and N were influenced by ICL and NE compared with corn–soybean without grazing or short-term oats/CC/grazing treatments. Integrated crop–livestock system and NE enhanced C and N concentrations in macroaggregates as well as in microaggregates. It is concluded that ICL and NE systems are sustainable prospects in enhancing overall soil health. Integrating crop and livestock improved physicochemical and microbial properties compared to the traditional corn–soybean system.


Results
Aggregate size distribution. Significant differences in ICL, NE, and CNT were observed among aggregate size distributions at sites 1, 2, and 3, however, at site 4, only 1-2 and 2-4 mm aggregates responded to the grazed pasture, oats, and CC treatments (Table 1). Integrated crop-livestock system at sites 1 and 2 resulted in a significantly (P ≤ 0.5) greater fractionation of > 4 mm aggregates than the NE and CNT. Moreover, ICL enhanced 0.053-0.25 mm aggregates at all three sites compared with NE and CNT. A higher proportion of 1-2 mm and 2-4 mm aggregates was observed under grazed pasture mix than the oats, oats + CC, and oats + CC + grazing treatments at site 4.
Integrated crop-livestock system and NE enhanced (P ≤ 0.05) MWD at site 1 by 28% and 29%, respectively, compared with CNT. Similarly, for site 2, average MWD was significantly greater for ICL and NE systems compared with CNT (Table 2). However, at site 3, ICL system showed 44% and 28% significantly higher MWD than NE and CNT, respectively. Natural ecosystem resulted in higher GMD than ICL and CNT at site 1, however, it did not differ significantly from ICL at site 2. Conversely, NE showed significantly lower GMD than ICL at site 3. The MWD and GMD did not vary statistically among any treatments at site 4. Integrated crop-livestock resulted in significantly higher (~ 47%) glomalin related soil protein than CNT at site 1 (Table 2). However, no statistical difference was observed among any treatments at the rest of the sites. www.nature.com/scientificreports/ Aggregate-associated carbon and nitrogen. Integrated crop-livestock system at site 1 showed a significantly higher SOC in > 4 mm aggregates with mean value of 57.1 g kg −1 aggregate, which was 36% and 49% greater than NE and CNT treatments, respectively (Fig. 1). For the aggregates 1-2, 0.5-1, 0.25-0.5 mm, ICL and NE showed significantly higher C than the CNT at site 1. Though at site 2, SOC did not vary significantly between ICL and NE treatments for aggregates > 4, 2-4, 1-2, 0.5-1, and 0.053-0.25 mm sizes, however, ICL had ~ 41% higher SOC than NE for 0.25-0.5 mm aggregates. At site 3, NE treatment performed considerably better than ICL and CNT for aggregates ranging between 0.25 to > 4 mm. Oats, CC, and grazing treatments at site 4 did not reveal any significant differences (P ≥ 0.05) in aggregate-associated SOC. At site 1, ICL system showed 33-44% higher N than NE and CNT for > 4 mm aggregates, whereas for aggregates 1-2, 0.5-1, 0.25-0.5, 0.053-0.25 mm, ICL and NE had significantly higher N than CNT (Fig. 2). At site 2, ICL had significantly higher N in > 4 mm and 0.25-0.5 mm aggregates than NE and CNT. Whereas at site 3, NE exhibited significantly higher N than ICL and CNT for all aggregate sizes. The only significant difference in soil N at site 4 was observed between oats and grazed pasture mix with mean values of 2.04 and 1.66 g kg −1 aggregate, respectively.
Overall, ICL and NE systems appeared more stable and favored the greater C and N associated to all aggregate size fractions, followed by CNT, where C and N were enhanced only in > 4 mm and 2-4 mm than the other aggregate sizes (Fig. 3a,b). On the contrary, cropping systems with oats and CC and grazing practices did not depict greater variations in C and N in different aggregate sizes. It was observed that aggregate-associated C and N are following similar trends under different cropping systems, thus also suggesting the tendencies to influence each other positively.
Soil microbial community composition. Natural ecosystem at site 1 showed significantly (P ≤ 0.05) greater total PLFA, total bacterial biomass, gram (+), gram (−) bacteria than CNT, and it did not vary significantly than ICL (Table 3). At site 2, NE showed significantly higher AM fungi and fungi/bacteria ratio than ICL and CNT. Natural ecosystem at sites 2 and 3 resulted in greater total PLFA than ICL, however, it did not vary significantly than CNT. Grazed pasture mix at site 4 depicted a total PLFA of 40.81 nmol g −1 soil, which was significantly greater compared to other treatments.
Principal component analysis. Principal component analysis was performed for measured soil physicochemical and soil biological properties i.e. aggregate size distribution, aggregate-associated SOC and N, glomalin, actinomycetes, gram (+), gram (−), AMF, total bacteria, total fungi, fungi/bacteria, and total PLFA (Fig. 4). www.nature.com/scientificreports/ www.nature.com/scientificreports/ The principal components (PC) 1 and 2 explained 33% and 22% of variation, accounting for the total variation of 55%. Different groups formed across the two principal component axes indicated the effect of different management systems on soil atributes. Integrated crop-livestock and NE were the systems with greater C and N fractions associated with different aggregate sizes and enhanced microbial biomass abundance. This represented the major influence of ICL and NE on soil health parameters compared to traditional corn-soybean without grazing or CC and short-term experiment with CC and grazing. Our results showed that soil aggregation and SOC and N in aggregates are related, therefore, may result in aggregate formation and protection of soil organic matter. Soil microbial composition, especially total PLFA, total bacterial and fungal biomass were greater under ICL and NE than CNT, grazed pasture mix, oat + cc, oats + CC + grazing, and oats, indicating that preservation of macroaggregates especially under long-term no-tilled or undisturbed systems positively influences the microbial population. Pearson's correlation illustrated the positive significant (P ≤ 0.05) influence of glomalin on MWD (r = 0.31) and GMD (r = 0.30); and AM fungi on SOC (r = 0.28) and GMD (r = 0.36) (Table 4). Overall, SOC had significant  www.nature.com/scientificreports/ (P ≤ 0.05) and positively correlated relationship with soil N and microbial population, e.g., AM fungi, gram (+), gram (−), actinomycetes, fungi/bacteria, total bacterial biomass and total PLFA. It is also noted that increased soil pH can negatively influence SOC and gram (−) bacteria; and higher soil bulk density (BD) can reduce gram (−) bacteria, actinomycetes, total bacterial biomass, and total PLFA. A strong negative correlation between BD and SOC (r = − 0.91) indicated that a reduction in soil C possibly due to increased decomposition can make soil more prone to compaction.   which is attributed to the presence of animals resulting in higher organic matter contents of total C and N fractions that can significantly enhance soil health over time 32 . Moreover, well-aggregated soils as observed under ICL (> 4 mm) at site 1 and NE (2-4 mm) at site 2 have a greater potential of retaining their structure and may have higher macropores, which facilitate sustained root growth than soils with low aggregation such as under CNT (corn-soybean without grazing or CC) in this study. It also explains the significance of ICL systems with no-tillage and undisturbed grassland, where the formation of stable macroaggregates may occur as a result of incorporation of plant residues, stimulation of root exudates and increased biological activity. Furthermore, it was noticed that ICL system not only enhanced the macroaggregates but accentuated the presence of microaggregates due to persistent binding agents, which are critical in SOC protection against microbial decomposition. When integrating grazing livestock into crop rotation, soil aggregation is typically improved under moderate and controlled grazing than the high intensity grazing systems 33 . Compared to the long-term sites (> 30 years), short-term site 4 did not result in discernible effects of grazing or CC on soil aggregation. However, within this short-term study, grazed pasture mix was able to enhance aggregation of 1-2 and 2-4 mm sized aggregates compared with oats, oats with CC, oats with CC and grazing. To observe the influence of CC and grazing on > 4 mm or < 1 mm aggregate sizes, a longer study period than 3 years might be required.

MWD, GMD, glomalin-related soil protein.
Overall, the observed higher MWD under ICL and NE than CNT showed the bonding of minerals and aggregates due to increased microbial activity and root exudation that improved overall aggregate stability 6 . Higher MWD for integrated systems than continuous corn monoculture has been reported previously 6 . The higher GMD at either ICL or NE system can be explained due to natural reconsolidation of soil than the CNT where soil can be disturbed due to conventional tillage 34 . Integrated crop-livestock at site 1 enhanced glomalin-related soil protein by 1.5 times when compared with CNT. However, this statistical difference (P ≤ 0.05) was only observed for 1 out of 3 long-term sites. Similar results were reported in previous study 35 , where glomalin content was not influenced under 15-19 years old grazed pastures at 0-20 cm. Although glomalin serves as an important aggregate binding agent and represents organicallybound soil N in agroecosystems, insignificant results suggest that independent of glomalin, accumulation of soil C can be more important in maintaining stable aggregates under long-term ICL systems.
Aggregate-associated C and N. Nearly constant C and N concentration was observed in aggregate size classes among soil under ICL and NE at sites 1 and 2, while relatively lower concentration was detected in macro and microaggregates under CNT. Large macroaggregates (> 4 mm) under ICL at site 1 resulted in 1.3-1.5 times significantly higher SOC concentration than NE and CNT. The greater concentration of SOC and N in ICL and NE is attributed to the lack of soil disturbance, crop residue retention, and rhizodeposition, which reduces macroaggregate turnover rate 14 . At site 3, NE enhanced aggregate-associated C and N concentrations and performed significantly better than both ICL and CNT treatments. The higher C and N accrual in the NE than ICL and CNT, especially at site 3, can be due to massive root systems, long-term establishment and absence of cultivation, which contributes to enhanced soil quality, while reducing nutrient vulnerability to loss by oxidation 18,36 . For short-term study at site 4, insignificant differences in aggregate-associated SOC suggested that longer study period of at least > 5 years is required for SOC to respond to grazing and cover crop management. The higher total N under ICL and NE can also be due to the presence of legumes, and brassicas in CC, which are effective at recycling N and may have helped in scavenging N. www.nature.com/scientificreports/ An overall increase in C and N cycling under ICL and NE systems has been attributed to ingested pasture being converted into urine and manure. Under these systems, livestock catalyze nutrient cycling by breakdown of complex plant molecules, greater soil incorporation and decomposition of plant residues and soil organic matter, which can maintain or even improve soil fertility by production of organic acids such as fulvic and humic acids 6,8,19 . Moreover, grazing stimulates the carbohydrate exudation from grass roots, which is mostly composed of polysaccharides, a C-O alkyl source 37 . The enhanced C concentration under ICL and NE can also be associated with higher MWD. Integrated system cool-season pasture and winter CC had significantly higher total C and N than the non-integrated continuous corn in previous study 6 . The results from another integrated system study 7 showed that soybean and oat-Italian ryegrass CC increased total C (1.16 Mg ha −1 yr −1 ) and N stocks (0.12 Mg ha −1 yr −1 ) under 7 year study period. It is previously reported that ICL system contains labile organic matter pools 10,38 , subsequently showing higher C stocks and greater root densities near soil-surface, which promotes aggregate-associated C stabilization 18,39,40 , higher infiltration rates, thus providing likely benefits to soil function linked to erosion control and soil water relations 41 .
Soil microbial community composition. Total bacterial biomass, AM fungi, and PLFA were enhanced under NE, which can be result of accumulation of organic residues and higher pasture root mass 7,32 , pasture being grazed can promote exudation of organic compounds by roots, serving as energy sources for microorganisms. The consistent increase in microbial population under NE can also be result of increased SOC and N, however, the same does not hold true for ICL system, where despite observing greater SOC and N, a significant decrease in the microbial population at site 2 was noticed. The enhanced total PLFA under NE system at site 2 is due to concomitant increase in AM fungi, gram (−), fungal/bacterial ratio, and total bacterial biomass compared to ICL. The fungal to bacterial ratio was reduced under ICL compared to NE at sites 1 and 2, pertaining to relatively low abundance of the fungal fatty acid 18:2ɷ6 in grazed system as compared to unmanaged grassland. This finding corroborates the notion that livestock-grazing systems contain bacterial-based decomposition channels and are mostly dominated by gram (+) bacteria and that the fungal population is comparatively more important in decomposer food-webs of native grasslands. These results coincide with previous studies 42,43 . Moreover, the increase in fungal to bacterial ratios under NE system in contrast to ICL at sites 1 and 2 can relate to modifications in soil health with C sequestration, as fungal populations incline towards higher C assimilation proficiencies and greater storage of metabolized C than bacterial populations 9,44 . The grazing intensity also plays a significant role in bacterial and fungal presence. It is previously reported that high grazing intensity had greater bacterial PLFA concentration than the low grazing counterparts in grassland systems 45 . It is considered that under heavily grazed sites in grasslands, bacteria-based energy channels of decomposition dominate other microbial communities, while fungi can successfully enable decomposition in both slightly grazed and non-grazed systems 43 . Grazed pasture mix at short-term study site 4 showed 12-21% higher total PLFA than the oats, oats with CC, oats with CC + grazing systems. It is also possible that this increased total PLFA at site 4 under grazed pasture mix contributed to enhancing the 1-2 and 2-4 mm sized aggregates compared to other treatments. It indicated that though physicochemical properties can take longer (> 8-15 years) in significantly responding to changes in management systems, soil microbial community and structure may show a rapid response (~ 3 years), thus it can be used as an early indicator while assessing the variations in soil health 18,46 .
Overall, NE exemplified the undisturbed grazed mixture with a greater microbial population at sites 1, 2, and 3, when compared to other agricultural systems. Our findings coincide with previous studies where pasture systems performed better than the agricultural soil, in terms of, showing greater microbial biomass and fatty acid signatures related to bacterial and fungal populations, which is mostly attributed to greater surface coverage and absence of tillage practices in pasture systems 9,47,48 . Lower soil microbial communities under ICL system than native Cerrado pasture have been found previously because of reduced soil porosity and macropore continuity resulting in restricted gas diffusion and water movement 18 .
Although the AM fungi abundance was not significant for sites 3 and 4, and significantly lower for ICL system than NE at sites 1 and 2, it should be taken into consideration that FAME analysis cannot reflect species-level changes for fungi and/or bacteria and the variations in microbial community structure for ICL system can be due to changes in abundance and distribution among microbial groups. For example, in a previous study 9 , while increased bacterial population was observed for continuous cotton compared to the ICL system, however, pyrosequencing for bacterial diversity assessment demonstrated disparities between both systems, where greater Proteobacteria was seen under ICL system than continuous cotton. Numerous factors such as degree of disturbance, pH level, bulk density, porosity, soil water content, C and N distribution, and residue positioning regulate the amount of bacterial and fungal biomass in agroecosystems 18,49 . Arbuscular mycorrhizal fungi are responsible for formation of macroaggregates (> 0.25 mm) by producing a glycoprotein called glomalin, which is present abundantly in natural and agricultural systems. However, increased grazing intensity, use of excess fertilizers and fungicides can directly or indirectly reduce mycorrhizal population by influencing soil organisms accountable for converting soil organic matter into plant nutrients 38 . Animals may also cause moderate soil compaction affecting the fungal biodiversity and soil pore space 6,38 . Relationship among measured soil properties. Based on PCA, it is derived that integrated croplivestock and natural ecosystem of native grassland can provide substrate for the microbial composition and enhance aggregate-associated C and N fractions. A positive correlation between SOC and microbial communities suggested the inclination of microbes to affect SOC and N turnover and vice-versa through interaction with crop-livestock grazing, vegetation, and soil properties. Fungi exhibited insignificant responses to changes in soil pH and bulk density than bacteria because chitinous cell walls make fungi more resistant and resilient to variations in soil conditions 50 www.nature.com/scientificreports/ conditions due to pH and increased bulk density, which has previously been observed in other studies 52,53 . A significant negative correlation between bulk density and SOC, N, gram (−), actinomycetes, total bacteria, and total PLFA indicated that the microorganisms influenced the soil compaction and related SOC and N. Moreover, positive correlation between SOC and microbial composition suggested that microbes can influence C sequestration in the soil via a shift in community structure. Microbial composition is influenced by soil C, whereas N is a critical biogenic element that improves microbial growth and their ability to utilize soil C 54 .

Materials and methods
Experimental site and treatments. The study was conducted in South Dakota at four locations consisting of three on-farm sites: site 1 (Salem, SD); site 2 (Bristol, SD); site 3 (Bristol, SD); and one experimental site 4 (Beresford, SD). The study sites have Köppen climate classification of warm humid continental weather with mean annual temperature of 6 °C and annual precipitation of 710 mm. Detailed site and treatment description is presented in Table 5. Sites 1, 2, and 3 consisted of similar long-term (> 30-year) experiments of: (i) integrated crop-livestock (ICL) system corn-soybean-oats (Avena sativa L.)-CC with cattle grazing, (ii) natural ecosystem (NE) as grazed native pasture, and (iii) corn-soybean without grazing or CC as control (CNT). Site 4 is an experimental site with study duration of 3-years consisting of treatments: (i) oats, (ii) oats with CC, (iii) oats with CC + grazing, and (iv) grazed pasture mix. Natural ecosystem was present at producers' locations for about > 50 years and part of this area was converted to croplands (ICL and CNT) by the producer. Natural ecosystem was used as a reference to compare any improvement that may occur in soil properties under ICL or CNT systems. The field size was 40 ha, with ICL and CNT treatments located within 50-m distance of each other, and natural ecosystem area was located 100-200 m away from the study sites. From November to March, a group of Aberdeen Angus cattle (Bos taurus) were used for grazing of corn and soybean residues, and CC every year 31 . Pastures were grazed during May and October. The experimental site 4 was set up as a randomized complete block design with four replications of plot size 18 m × 36 m. From this site, four soil samples were collected from each plot and were composited for soil analyses. At producers' sites, three pseudo-replicates were established depending on the landscape position and soil characteristics. Four soil samples were taken from each established pseudo-replicate. Cover crop treatment consists of legumes, grasses, and brassicas such as pea, sorghum, radish, lentil, cowpea, turnip, and oats at four sites. The study complies with local and national regulations. All methods were carried out in accordance with relevant guidelines and regulations and all experimental protocols were approved by the Institutional Animal Care and Use Committee of South Dakota State University. Seeds used in the study were supplied by Millborn Seeds Inc (Brookings, SD) and are property of South Dakota State University.
Soil sampling and measurements. Aggregate size distribution. Soil samples were collected during summer 2019 at 0-5 cm depth using a push probe auger. The wet sieving method of aggregate separation was performed as suggested by other studies 55 , and 56 . Air-died soil (100 g) was submerged in wet sieving vertical shaker apparatus (CSC sieve shaker, Fairfax, VA) adjusted to 90 rpm and 1.3 cm stroke length for 5 min on top of a set of 4, 2, 1, 0.5, 0.25-, and 0.053-mm sieves. Aggregate size distribution was determined by soil remained on each sieve and further dried at 105 °C oven temperature, weighed and expressed as a fraction of the initial mass of samples used. The mean weight diameter (MWD) was calculated using the following equation 57 : www.nature.com/scientificreports/ where xi is the mean diameter of aggregate size on each sieve (mm), wi is aggregate mass retained on each sieve (g) as a fraction of total dry weight of soil (g), and n is the number of aggregate size fractions. The geometric mean diameter (GMD) was calculated by the equation 58,59 : where xi is the mean diameter of each size fraction (mm) and wi is the fraction of total soil sample weight present in the size fraction i.
Aggregate-associated C and N. Aggregate size fraction retained on each sieve was ground and a 0.250-g soil sample was weighed in a tin foil cup for analyzing SOC and N in aggregates with a LECO TruSpec Analyzer (LECO Corporation, St. Joseph, MI) using the dry combustion method.
Phospholipid fatty acid (PLFA) analysis. Extraction. To assess microbial community structure, total soil lipids were extracted by shaking approx. Glomalin-related soil protein. Soil protein was extracted using the method provided by 60 with some modifications by 61 . Air-dried soil (3 g) was weighed into a glass tube and 24 ml of 20 mM sodium citrate buffer was added, followed by shaking at 180 rpm for 5 min. The tubes containing samples were autoclaved at 121 °C, 15 psi for 30 min and cooled. Two-ml of the slurry was retracted to a microcentrifuge tube and soil particles were separated by centrifuging the slurry at 10,000 × gravity. Two-ml of Thermo Pierce™ bicinchoninic acid (BCA) protein assay solution was added to 0.1 ml of clarified extract and was incubated for 1 h at 60 °C. Color development was read using BioTek spectrophotometric plate reader. Extractable soil protein content was determined by multiplying the extract's protein concentration by the extractant's volume and dividing by soil (g) used.
Statistical analyses. The SAS 9.4 version (SAS Institute, Cary, North Carolina, USA) was used with oneway ANOVA procedure to analyze the influence of different crop management systems on soil physical and biochemical parameters. Each site was analyzed separately where ICL, NE, and CNT treatments were considered as fixed effects. Fisher's protected least significant difference (LSD) was used for mean separation of treatments at a significance level of 0.05. Shapiro-Wilk's test was used to analyze the normality of the data and Levene's test was used to assess the homogeneity of variances. Principal component analysis (PCA) of the full dataset was analyzed using the PRINCOMP procedure of SAS 9.4. Pearson's correlation was performed using different soil physicochemical properties and soil microbial structure using PROC CORR statement of SAS version 9.4. The PCA was used to subgroup experimental treatments based on the measured soil parameters by generating the eigenvectors of these parameters called loadings and principal component (PC) scores for each unit. Each eigenvector loading indicates the direction and magnitude of association between soil properties and treatments 62 . www.nature.com/scientificreports/

Conclusion
This study evaluated the influence of integrated crop-livestock management, cover crops, and grazing compared to the traditional corn-soybean rotation on selected soil health indicators focusing mainly on aggregate size distribution, MWD, GMD, glomalin-related soil protein, aggregate-associated C and N, and soil microbial community structure. The performance of these measured soil parameters under integrated crop-livestock system was also compared with > 50 years old native grassland used as a reference. The overall significantly greater performance was observed in the soil parameters at sites 1, 2, and 3, indicating that the long-term establishment of producers' sites for > 30 years was efficient and showed considerable differences in soil health contrasting to the experimental site which had been established for a comparatively shorter duration of 3 years. In particular, ICL resulted in improved aggregate fraction with greater > 4 mm and 0.053-0.25 mm aggregates compared with NE and CNT, attributed to animal presence, incorporated manure, and improved microbial diversity. Insignificant differences in glomalin protein at sites 2, 3, and 4 among all treatments represented that independent of glomalin, soil C could play an important role in binding aggregates. Improved C and N cycling was observed not only in macroaggregates but also in microaggregates under ICL at site 1 and both ICL and NE at site 2 when compared with corn-soybean control. Natural ecosystem performed well in terms of increasing total PLFA at all three producers' sites, which can be due to higher vegetation cover, higher pasture root mass, absence of cultivation, hence higher root-associated organic compounds enhanced due to grazing. Phospholipid fatty acid analysis showed that ICL and NE at site 1 had higher total PLFA, total bacterial biomass, gram (+), gram (−) bacteria than corn-soybean rotation, suggesting the improvement in the soil microbial structure by the addition of livestock into the crop production system. The response of grazed pasture mix under short-term experiment (~ 3 years) was inconspicuous in terms of variations in aggregate-associated C and N, however, it did increase the 1-2 and 2-4 mm aggregate distribution compared with the oats/CC/grazing practices, possibly due to observed higher total PLFA. The observed significant differences in aggregation and microbial composition at the experimental site 4 suggest that these may serve as physical and biological indicators for short-term experiments and can better detect changes in soil than C and N variations as a result of different management systems. Integrated crop-livestock system may have similar or greater benefits as a long-term natural ecosystem in terms of increasing macroaggregate and microaggregate fraction along with associated C and N cycling, total PLFA, total bacterial biomass, gram (+), gram (−) bacteria compared to traditional corn-soybean without grazing or cover crops.